Entanglement verification with realistic measurement devices via squashing operations 
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Many protocols and experiments in quantum information science are described in terms of simple 
measurements on qubits. However, in a real implementation, the exact description is more difficult, 
and more complicated observables are used. The question arises whether a claim of entanglement 
in the simplified description still holds, if the difference between the realistic and simplified models 
is taken into account. We show that a positive entanglement statement remains valid if a certain 
positive linear map connecting the two descriptions — a so-called squashing operation — exists; then 
lower bounds on the amount of entanglement are also possible. We apply our results to polarization 
measurements of photons using only threshold detectors, and derive procedures under which multi- 
photon events can be neglected. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.65.Ta 
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I. INTRODUCTION 



According to Asher Peres, entanglement is "a trick that 
quantum magicians use to produce phenomena that can- 
not be imitated by classical magicians" Because of 
the key role of entanglement in applications lots of effort 
is put into realizing this fragile resource in the lab, for 
example via parametric down-conversion (PDC) sources 
or with ion traps, to only name a few. In a real experi- 
ment it is of course desirable to unambiguously verify the 
creation of entanglement, and in fact many different op- 
erational tools have been developed over the past years to 
achieve this task, cf. Ref. [|J for a review. A reliable en- 
tanglement verification has to satisfy a few crucial criteria 
3 ; most importantly the verification method should not 
rely on assumptions from the entanglement generation 
process, but instead on the information acquired about 
the system via measurements. Moreover the obtained 
data should be considered under a worst case scenario. 
That is, in the spirit of Ref. Q, the test is only consid- 
ered to be affirmative if, in the limiting case of an infinite 
number of experimental runs, the data exclude compat- 
ibility with all separable states. This viewpoint is even 
essential for certain tasks like quantum cryptography Q . 

In any case, it is typical to allow one basic ingredi- 
ent: since one usually believes in quantum mechanics, it 
is common to assume that one is equipped with an ac- 
curate quantum mechanical description of the employed 
measurement devices; the actual testing or the (experi- 
mental) characterization of a measurement device is any- 
way often combined with other assumptions concerning 
the generated state [1, 0, 9 ■ Note that if one does not 



restrict oneself to this model then one can still use Bell 
inequalities for the verification. This leads to the known 
drawback that the entanglement of certain states can 
never be verified 0] and there is even the conjecture that 
complete classes of interesting entangled states may be 
undetectable [l(| . However this will not concern us here, 
and we always assume to have an operator set associated 
with the observed data, which allows us connect the data 
to quantum mechanical quantities. 

An example of a straightforward and hence quite often 
applied entanglement verification method, e.g., Ref. [Til ] . 
is the following procedure which we call the tomography 
entanglement test in the following: Since the useful en- 
tanglement might be confined to a low dimensional sub- 
space, e.g., the single photon-pair subspace of a PDC 
source or two very long-lived energy levels of two ions in 
a trap, one just performs a few different measurements to 
obtain tomography on this subspace. After several runs 
of the experiment one has collected enough data to recon- 
struct the underlying density operator on this subspace 
via some reconstruction technique. Note that here one 
employs the knowledge of the measurement description. 
In order to check for entanglement one just investigates 
whether this reconstructed density operator describes an 
entangled state or not. 

However, does one really verify entanglement via this 
method? The problem lies within the measurement de- 
scription, because such ideal measurements, as the ones 
used in the reconstruction mechanism, might not have 
actually been performed in the experiment. A good ex- 
ample is represented by the polarization measurement 
with two threshold detectors (see also Fig. [5]) , which is 
typically employed in photonic experiments. Apart from 
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usually acting onto several input modes at once, this de- 
vice does not even respond solely to the single photon 
subspacc, since such detectors cannot resolve the num- 
ber of photons. Hence the question arises whether one 
still verifies the entanglement if a more realistic measure- 
ment description is employed. It is the main purpose of 
this paper to study this question. Note that the afore- 
mentioned scenario often occurs, not because one is not 
aware of the more realistic model, but because an over- 
simplified measurement description is employed in order 
to ease the task of entanglement verification. 

Specific instances of the problems considered here have 
been investigated in several works in the literature. In 
Ref. 12j inequalities for the detection of entanglement for 
two qubits have been proposed, where the measurement's 
devices can be misaligned to a certain degree. Bell-type 
inequalities which are independent of the spectrum of the 
measured observablcs have been recently introduced in 
Ref. [l3l ]. Moreover, for an experiment with photons from 
atomic ensembles, an entanglement verification scheme 
which takes multi-photon events into account has been 
introduced [l4| an d implemented [l5j . 

In this paper we proceed along the following lines: In 
Section II we provide an example of a tomography entan- 
glement test which indeed leads to the wrong conclusion 
about the presence of entanglement under a small, phys- 
ical change of the employed measurement description. 

In Section III, we start to investigate under which con- 
ditions such mistakes can safely be excluded. In short, 
the entanglement verification process remains stable as 
soon as the considered set of operators are connected by 
a positive but not necessarily completely positive map, 
the so-called squashing operation. Similar relations be- 
tween different measurement schemes have recently been 
introduced in the context of QKD, cf. Ref. 0, an( i 
even other known verification methods can be cast into 
this framework. However, complete positivity of the con- 
nection map was required there. 

In Section IV we reformulate the existence of such a 
positive map into a necessary and sufficient condition 
which provides a particular intuitive solution for the to- 
mography entanglement test: The map exists if and only 
if each classical outcome pattern from the refined set of 
observables remains compatible with the oversimplified 
set of observables. 

Then, in Section V we prove that the aforementioned 
polarization measurement with threshold detectors along 
all three different polarization axes represents an example 
which indeed can only be linked to its single photon re- 
alization by a positive but not completely positive map. 
This analysis concludes that the tomography entangle- 
ment test which is typically employed for a PDC source 
[T§ | or even in multipartite photonic experiments fl9j (us- 
ing the single photon assumption) can indeed be made 
error-free if the (local) double click events are taken into 
account. 

In Section VI we consider the issue of entanglement 
quantification, proving that one can in principle get lower 



bounds on the entanglement of the physical state in terms 
of the entanglement of the operator that results from the 
local mapping between the observables. 

Finally, we conclude and provide an outlook on possi- 
ble further directions. 



II. AN EXAMPLE FOR ION TRAP 
EXPERIMENTS 

Let us first mention a simple, yet practically relevant 
example, which shows that the tomography entanglement 
test indeed can lead to a false conclusion about the pres- 
ence of entanglement if the structure of the observables 
is not properly taken into account. 

For a single 40 Ca-ion in a trap one typically consid- 
ers only the lowest two energy levels given by a lower 
level 1 5*) = |1) and the upper level \D) — |0) and treats 
them as the qubit (20j . Resonance fluorescence provides 
a mechanism to read out the occupation number of the 
energy levels: An electron in the \S) state is coupled to a 
higher energy level |P), and observing photons from the 
IS*) <-> \P) transition signals that the qubit was in the 
state \S). This overall process corresponds to a projec- 
tion onto the lower energy state and consequently allows 
to measure the a z Pauli, while the measurement along 
different directions is achieved by a local basis rotations 
prior to the a z measurement, cf. Ref. poj . 

In order to avoid too many measurements it is com- 
mon to measure the occupation probability only for the 
state \S), simply because for qubits the other probabil- 
ity equals p(D) = 1 — p(S) due to the normalization, and 
similar for the remaining basis settings. Suppose that one 
uses this measurement procedure to obtain tomography 
in order to verify the creation of entanglement between 
two separated ions in the trap. Consider now the example 
that the observed expectation values, abstractly denoted 
as Eij (p) and characterized by a noise parameter p {52| , 
may allow the reconstruction of the state 



p(p) = (l-p)\iP + )(iP + \+p~, 



(1) 



which is, by virtue of the PPT criterion, entangled for 
p < 2/3. 

However in practice the situation is more complicated 
since the ion is not a simple two-level system. To model 
this, one can add another energy level to only one of 
the ions, thereby enlarging the two-qubit system to a 
qubit-qutrit one. Without any additional information 
about the occupation number of this extra level, it is 
clear that the assignment p{D) = I — p(S) is not cor- 
rect any more. Consequently the observed data Eij (p) 
can only verify entanglement for the case p < 0.63. This 
can be checked by using the tools from Ref. [2l[ , in which 
the search for an appropriate separable state was phrased 
into a semidefinite program. Hence we have the inter- 
val p £ [0.63,2/3), for which the performed tomogra- 
phy entanglement test indicates the presence of entan- 
glement although with the more realistic model it does 



not. Though this region might be small this error can 
become important in the multipartite scenario, where 
current experiments just operate at the border of gen- 
uine multipartite entanglement [12, [H, [24[ ■ Concerning 
the experimental consequences, however, two facts are 
important: 

1. For experiments with ion traps it is known that the 
occupation probability for levels apart from the two 
logical states is very small, typically it is around 
10 -3 [53| . Given this additional measurement data, 
it is possible to provide a quantitative estimate of 
the resulting error in the used entanglement verifi- 
cation scheme, e.g., the mean value of an entangle- 
ment witness. For typical entanglement witnesses 
employed in those scenarios this error is far be- 
low the unavoidable statistical uncertainty, which 
is caused by the finite number of copies of a state 
available in any experiment. 

2. Note that the probabilities p(S) and p(D) of each 
energy level can be measured independently by ad- 
ditional local rotations, hence at the expense of 
more measurements. Then the resulting proba- 
bilities correspond to the unnormalized two-level 
state p rc d that is obtained from our modeled three- 
level system p tot by a local projection, i.e., p Te d = 
nptotll, with n = \S)(S\ + \D)(D\. As long as 
we prove entanglement of the two-qubit system 

= 1 ® Ilp^t 1 ® n, this also implies entangle- 
ment for the total state p^, since the projection 
is local. 

For instance, if one measures a witness like 
W = |00)(00| + |11)(11| - \x+x + )(x+x+\ - 

\x~x~){x~x~\ + \y + y + )(y + y + \ + \y~y~)(y~y~\, 

with = (|0) ± \l))/y/2 and |y±) = (|0) ± 

i|l})/\/2 [H, the mean value of W is nothing more 
than a linear combination of certain probabili- 
ties on the qubit space, and if the mean value is 
negative, the state p^ and hence p^ is entan- 
gled. This shows that additional dimensions of the 
Hilbcrt space alone do not invalidate the conclusion 
that the state is entangled when the measurement 
devices are characterized properly. 



III. POSITIVE SQUASHING OPERATIONS 

We are ready to formulate the problem that we solve 
throughout the subsequent sections. For each local mea- 
surement setup one has two different sets of ordered ob- 
servables; a set of simple target observables labeled as Tj 
with i = 1, . . . , n acting on the Hilbert space H.t which 
are used for the entanglement verification process or in 
the reconstruction mechanism, and a different set of so- 
called full operators denoted as Fi with i = 1, . . . ,n onto 
the Hilbcrt space TIf which represent the more realistic 
model of the actual observables in the experiment. In the 



3 



Full observable set {Fi} 




FIG. 1: Idea of the positive squashing operation: The given 
fall observable set {Fi} can be decomposed into a positive 
squashing operator A followed by a particular target observ- 
able set {Ti} such that the same expectation values Ei are 
obtained for all possible input states pF- 



above ion-trap example we considered the case of qubit 
target observables, while our full operators were acting 
on a qutrit system. 

Consider the case where in an experiment one mea- 
sures the expectation values of the full operators Fi but 
instead one interprets them as the expectation values of 
the target observables Tj. The question arises, whether 
this may lead to a false entanglement verification. In the 
following we provide a simple condition on the two oper- 
ator sets alone that excludes such a possibility, and hence 
guarantees the presence of entanglement. 

Suppose that both sets of observables are connected by 
a positive (but not necessarily completely positive) linear 
map A : C(Hf) — * C(TIt) which satisfies the following: 
the expectation value of any observable Fi with respect to 
an arbitrary input state pp is the same as the expectation 
value of the corresponding target operator Ti with respect 
to the output state of the corresponding map px = A(pp) 
(see Fig. [T]). That is, 

tr(p F F i )=tr[A(p F )T i \ (2) 

holds for any input state pv and for all? = 1, . . . , n. Using 
the adjoint map A^ : C(Ht) — * C(Hf) this condition can 
be rephrased into 

A^(T i )=F i (3) 

for all i = 1, ...,n, while the positivity requirement 
transfers also to the adjoint map A^ . 

Such a described connection between two different ob- 
servables sets is an extension of the notion of a squashing 
operation in the QKD context [lj| [ljj which differs from 
the present definition only by the extra condition of be- 
ing completely positive and trace-preserving. Because of 
these similarities we use the term positive squashing oper- 
ation in order to refer to map A (or its adjoint A^). Note 
that typically we consider the case of a trace-preserving 
map A (or unital map A^) such that density operators are 
mapped to properly normalized density operator; how- 
ever this requirement is not mandatory. An example of 
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a non-trace preserving, but positive map between opera- 
tor sets is given by the matrix of moments 54 1 (2f|; the 
only difference is that one must be careful with entangle- 
ment criteria on the target space that explicitly employ 
the normalization of the density operators (e.g., the com- 
putable cross norm or realignment criterion) , but one can 
also deal with this (2f| . 

The advantage of such a positive squashing operation is 
that the structure of separable states (from the full to the 
target Hilbert space) remains invariant, and hence any 
successful entanglement verification on the target space 
directly translates to a positive verification statement on 
the full Hilbert space: 

Proposition III.l (Entanglement verification). Let us 
assume that the two sets of local observables on Alice's 
side, labeled as {T 4 A } and {F^} respectively, are con- 
nected by a positive (not necessarily completely positive) 
unital linear map satisfying Eq. f3|) ; and a similar re- 
lation holds for Bob 's side. If the observed data verify en- 
tanglement with respect to the target observables T^(&T^ , 
then this data also proves the presence of entanglement 
for the full operator set <X> F^ . An analogous state- 
ment holds for more than two particles. 



Proof. For the observed data Eij one has the identity 
Ekj = tr^AB-F^ <8> if) = tr[A AB (/3AB)T i A ® Tf] due to 
the property of the squashing operation. For any sep- 
arable state on the full bipartite Hilbert space p^B — 
Y^kPkPA ® Pbi onc obtains 



sep * / se P\ 
7 AB : = A Ab(Pab) 



b), 



(4) 



which represents a valid (normalized) separable density 
operator on the bipartite target Hilbert space because 
of positivity of the corresponding (unital) maps, and is 
compatible with the observed data. Consequently, if one 
proves the incompatibility of the mean values of the Ti 
with all separable states on the target space, the density 
matrix on the full space must be entangled. Note that 
here one just needs positivity of Aa and Ab and not 
complete positivity. □ 

Note that a local squashing operation between opera- 
tor sets does not represent the most general map between 
bipartite observable sets that preserve the structure of 
separable states; however we neglect other options on 
behalf of the "locality" of this connection. Furthermore 
note that since we do not require for a completely positive 
map, it can happen that one obtains an unphysical (not 
positive semidefinite) density matrix on the target space; 
such an operator is then also incompatible with a sepa- 
rable state. However this situation can only occur for an 
entangled state on the full bipartite Hilbert space, hence 
the conclusion of the entanglement verification process 
remains unaffected. 

Finally, let us add that the precise state reconstruction 
technique needed for the tomography entanglement test, 



cither direct inversion of Born's rule or maximum likeli- 
hood estimation (2|| (although there are even problems 
associated with them [231), does n °t conflict with a pos- 
itive but not completely positive squashing operation. If 
the corresponding operator on the target space is posi- 
tive semidefinite both reconstruction techniques deliver 
the same operator (in the limit where one really obtains 
exact knowledge of the expectation value). Because any 
separable state is represented by a valid separable target 
state this excludes the possibility that a separable state 
is mapped to an entangled state by the reconstruction 
process. In the case of an unphysical "entangled" target 
operator a direct inversion of Born's rule one would di- 
rectly "witness" the entanglement (HJ. In contrast the 
maximum likelihood method produces the closest posi- 
tive semidefinite operator (2?| (with respect to the likeli- 
hood "distance"), hence an unphysical, entangled target 
state can be mapped to a separable state via this re- 
construction technique and thus escapes the tomography 
entanglement test. But this does not bother us here, be- 
cause some entangled states are missed anyway due to 
the simplified operator set. 



IV. CRITERIA FOR THE EXISTENCE OF A 
POSITIVE SQUASHING OPERATION 

In this section we investigate which requirements need 
to be fulfilled by the two different operator sets in or- 
der to be connected by a positive squashing operation. 
There are, of course, different ways how one can tackle 
this problem: One method, in close analogy to that of 
Ref. I17f. is to employ the Choi-Jamiolkowski isomor- 
phism |28l. l29l . |30| between linear maps and linear opera- 
tors. This isomorphism transforms positive maps rcdinto 
entanglement witnesses, or more precisely into linear op- 
erators that are positive on product states, while the re- 
quirements from Eq. ([3]) change into a set of linear equa- 
tions that constrain the allowed form of the entangle- 
ment witness. For an explicit solution to this reformu- 
lated problem one first solves these linear equations and 
afterwards tries to choose the remaining, undetermined 
parameters of the operator in such a way that it meets 
the entanglement witness condition. 

However, we take a different path that provides a clear 
interpretation for the existence of such a positive linear 
map and which is also employed in the later part to prove 
the positive squashing property for the polarization mea- 
surements. 

Equation (JSJl directly allows us to read off a necessary 
condition: it states that for each physical density op- 
erator p-p in the full Hilbert space there exists a valid 
density operator A(/?p) (if A is trace-preserving) in the 
target space such that both operators assign the same 
expectation values for the considered observables. Hence 
all possible expectation values Ei that can in principle be 
observed on the full Hilbert space must remain physical 
with respect to the target observables. As we will see, 
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this condition becomes also sufficient if the target opera- 
tors Tj with i = 1, . . . , n provide a complete tomographic 
set. Thus, in combination with Prop. IIII.ll we have the 
following solution for the question posed in the introduc- 
tion: The tomography entanglement test is error-free as 
long as the full local observables on Alice and Bob's side 
can only produce measurement results which are also con- 
sistent with the local target, or reconstruction observables. 

For the following proposition we need to define the set 
of possible physical expectation values associated with a 
given set of observables, defined as 

S F = [-E G M n |there is a p G T>(H F ) such that 

Ei = tr(pFi) for all i = 1, . . .re}, (5) 

and a similar definition for the operator set on the target 
system <Sp- Concluding we have the following character- 
ization: 

Proposition IV. 1 (Existence). The set of full observ- 
ables {Fi} and the tomographically complete set of target 
observables {Tj} are related by a unital squashing opera- 
tion if and only if it holds that Sf Q St ■ 

Proof. One direction of the proof is clear: Suppose that 
there exists a positive trace-preserving squashing opera- 
tion A. For any E G Sf we must have a density operator 
Pf such that one obtains Ei = tr(pp-Fi) = tr[A(pp)Tj]. 
Because of the properties of the corresponding map we 
receive a valid target density operator pp := A(/Op) which 
provides the same expectation values E, hence E G St- 
This concludes the first direction of Sf Q St- 

For the other direction, we employ the fact that the set 
of target operators are tomographically complete and the 
set inclusion Sf Q St to explicitly write down the posi- 
tive squashing operation. First note that for a given set 
of physical expectation values E G St, the corresponding 
target density operator is uniquely determined by a di- 
rect inversion of Born's rule, 72-t : E i— > pr(E), i.e., by a 
linear reconstruction mechanism that maps the expecta- 
tion values to its explicit density operator. Moreover for 
a given full density operator pp the corresponding expec- 
tation values are already determined, which is described 
by the linear map M-p : pp i— ► E. Combining these two 
maps according to 

A = K t °Mf (6) 

provides the squashing operation: That is, for a given 
input state pp one first computes the expectation values 
Ei of the full operator set and then uses these values in 
the reconstruction algorithm (that depends on the tar- 
get operators) to obtain the corresponding target output 
state. The set inclusion guarantees that any valid full 
density operator is mapped to a valid target state, hence 
the described map is already positive. Since both maps 
in the decomposition are linear the overall map is linear 
as well. □ 



In a concrete example the proposition of course only 
helps if one obtains knowledge on the sets iS; although 
this is by far not a trivial task, one can employ approxi- 
mation techniques for a special set of observables or even 
a hyperplane characterization for the exact determina- 
tion, see Ref. [3l| for more details. 

If the set of considered observables on the target space 
is not tomographically complete, then, in order to estab- 
lish the existence of a positive mapping between the two 
sets of operators, one can still invoke Proposition IIV.1I 
with some caution. Indeed, one has to check whether it 
is possible to extend the two sets by some choice of addi- 
tional target and full operators, so that the target set is 
tomographically and the two sets of physical expectation 
values — which depend on the choice of the extensions — 
satisfy the condition of Proposition IIV.1I 

Finally, let us note that one can also characterize a 
completely positive map via such a set inclusion require- 
ment if one adds an additional reference system R of 
dimension equal to that of the full space (or of the target 
space, in the case the dual map) on each side, because 
complete positivity of A just means that idpt ® A is posi- 
tive for such a reference R. For an actual investigation of 
such a completely positive map, however, the formalism 
of Ref. [ijj seems more appropriate to us. 



V. EXAMPLE: POLARIZATION 
MEASUREMENTS 

In this section, we apply the developed formalism to a 
physical relevant measurement setup. We draw our at- 
tention to polarization measurements onto a two-mode 
system by using only threshold detectors, i.e., such de- 
tectors cannot resolve the number of photons. More 
precisely, as shown in Fig. [21 the incoming light field 
is separated according to a chosen polarization basis 
/3 G {x, y, z} via a polarizing beam splitter, followed by a 
photon number measurement on each of those modes by 
a simple threshold detector. Hence in total four differ- 
ent outcomes can be distinguished: no click at all, only 
one of the detector clicks or both of them trigger a signal 
and produce a double click. Because of its great sim- 
plicity this measurement device appears quite frequently 
in quantum optical experiments which employ the polar- 
ization degree of freedom (for an overview see Ref. [32|V 
It turns out that this measurement device provides, if 
measured along all three different basis settings, a non- 
trivial example for the difference between a positive and 
a completely positive squasher. This means that the cor- 
responding map A can only chosen to be positive but not 
completely positive. 

Next let us specify which observable sets should be 
connected by the squashing operation; see also Ref. [l~7l ]. 
For each chosen polarization basis (3 the different mode 
measurements are denoted as M^p with the label i G 
{vac, 0,1, d} for all classical outcome possibilities: vac- 
uum, click in "0", click in "1", and a double-click. The 
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FIG. 2: Schematic setup of the considered polarization mea- 
surement: Via quarter (QWP) and half wave plate (HWP) 
one can effectively adjust the polarization basis (3 of the cor- 
responding polarizing beam splitter (PBS) according to the 
basis {±45°, O / (3, H/V} that we label as {x, y, z}. 



target measurements are chosen such that they justify 
a single photon description: each click event is inter- 
preted as a single photon state, hence as target mea- 
surements one selects the same measurement device, only 
that it just acts on the single photon subspace and the 
vacuum component. In order to achieve the squashing 
property, the double click events must be taken into ac- 
count, since such events are clearly incompatible in a 
(perfect) single photon interpretation, but they never- 
theless contribute to the normalization. One can incor- 
porate this effect by a particular post-processing scheme 
that represents a sort of penalty for double click events. 
A common scheme, originally introduced for the QKD 
context in Ref . [33|, [34| , consists of randomly assigning 
each double click event to one of the single click out- 
comes. This particular set of processed measurement op- 
erators becomes the exact set of full operators {F iy p} with 
i e {vac, 0, 1} and (3 G {x, y, z} and with the relation 
Fiji = Mi j 3 + l/2M d ,0 with i = 0, 1 for all [3. 

Let us start with a perfect single-polarization mode 
description of the full operators; imperfections like fi- 
nite effici ency or dark counts are considered later (see 
also Ref. j33[). The "no click" outcome is independent 
of the chosen polarization basis and becomes F vac / 3 = 
|0, 0) (0, 0| . All other observables are block-diagonal with 
respect to the photon number subspace, i.e., F^p = 
X)n=i ^Tp an d for a fixed number of photons we have 



!„ + (-!)* |n,0)fl(n,0|-|0,n) /J <0,»| , (7) 



with i = 0,1. Here \k,l)g denotes the corresponding 
two-mode Fock state in the chosen polarization basis f3 
(e.g., for n = 3 the state |2, 1) = \2H, IV) describes 
a system with two horizontally and one vertically po- 
larized photon) and l n represents the identity operator 
onto the n-photon subspace, which appears because of 
the chosen post-processing scheme. This perfect single- 
polarization mode description is also employed for the 



target operators, however only acting on the vacuum 
7vac,/3 = 1 0, 0} (0, 1 or on the single photon subspace 

T *,/3 = K,/3 with * = 0,1. 

Let us further comment on these observable sets: Note 
that if one selects the following standard basis for the 
single photon subspace |1,0) 2 = |0) and |0, 1) 2 = |1), 
then each difference of the single click outcomes equals to 
a familiar Pauli operator, i.e., up = Tq ^ — T\ p for all (3. 
Hence each of the single click operators T^p with i = 0, 1 
corresponds to a projection onto one of the two different 
cigenstates of the related Pauli operator a p. Furthermore 
the corresponding difference between the full observables 
Fp = Fqjj — i<i )( g is again block-diagonal and each n- 
photon part is given by 



^ = Fqp - F™p = \n, 0} p (n,0\ - \0,n) p(0,n\, 



(8) 



according to Eq. 0. Note that these observables are 
also accessible with a different polarization measurement 
that only uses a single threshold detector (5(| and which 
has alternativel y b een employed for polarization experi- 
ments, cf. Ref. [Ill ] - 

The following theorem proves the positive squashing 
property between the two given sets of observables; how- 
ever it also applies to the other measurement description 
of Ref. [ll|. 

Theorem V.l. There exists a positive, but not com- 
pletely positive unital squashing operation for the 
operator sets {Ti t p} and {Fi t p}, i.e., hJiT^p) = F^p. 
Therefore, the interpretation of the {Fi t p} as single pho- 
ton measurements {Ti t p} does not invalidate the entan- 
glement verification scheme. 

Proof. First let us point out that the existence of a com- 
pletely positive squashing operation has already been 
ruled out in Ref. [17j . 

In order to prove the existence of a positive squashing 
operation we only need to focus on the "click" events, 
since the vacuum part can be directly removed by a pro- 
jection discriminating between the vacuum and all other 
Fock-states. Note that it is sufficient to prove the squash- 
ing operation for a complete set of linear independent tar- 
get operators only, because other linear dependencies are 
implicitly present in the linear map. In short, it is equiv- 
alent to prove a unital squashing operation A\o~p) = Fp 
for all (3 £ {x,y,z}, where Fp is the described difference 
between the click outcomes of the full observables. 

Since we only concentrate on the single photon sub- 
space we are equipped with a full tomographic set and 
hence can readily apply Prop. IIV.ll such that it remains 
to prove Sf C St- Since each full observable is photon 
number diagonal one obtains that Sf is given by the con- 
vex hull of all n-photon sets Sp, i.e., the set of physical 
expectation values on an n-photon state. Hence we need 
to verify that each n-photon state can only produce ex- 
pectation values which are also compatible with a single 
photon state, i.e., Sp C Sp = St for all n > 1. The set 
Sp directly equals to the familiar Bloch sphere. Hence 
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we prove existence of a positive squashing operation if we 
can show that 



]T [tr(p^")] 2 <l 



(9) 



holds for all n photon density operators p, and for all 
photon numbers n > 1. 

In order to simplify the analysis in the following, each 
operator Fjf can be regarded as an operator acting onto 
an n-qubit space. Indeed, the n-photon Hilbert space 
Hp = <C n+1 is isomorphic to the symmetric subspace 
Sym(H„) of an n-qubit system H n = ((D 2 )®". Using the 
given standard basis definition one obtains for example 



F? = 



-|1>(1| 



50 ri 



(10) 



while for any other operator FJj one just replaces the 
states |0), |1) with the eigenvectors of the corresponding 
Pauli matrix erg. 

Expanding these operators in a multi-qubit basis de- 
livers 



17" 



1 + cr 



E EX* 



j®(n-j) 



(11) 



j odd 7T 



in which denotes the sum over all possible permuta- 
tions 7r(-) of the subsystems that yield different terms. 

Next, we exploit the result from Ref. [35[ that for odd 
j every quantum state p. hence also each state on the 
symmetric space, satisfies 



j3=x,y,z 



with the abbreviation 



Tr(af)) =tr 



piT\(J 



i(n-j) 



(12) 



(13) 



This inequality is based on the property that the ob- 
scrvables 7r(er| i ® 1®(»-J)) with p g {x,y,z} have all 
eigenvalues equal to ±1 and anti-commute pairwise [57| . 
Note that this identity holds for all occurring j and for 
all possible permutations it. Consequently one obtains 



E Mpf ( 



n 2 



2 2r, 



i E E 

odd tt-tt' 



(14) 



<1. 



where the inequality (together with the Cauchy-Schwarz 
inequality) was used to upper bound each term in the 
squared bracket by 1. For the final result one needs to 
count the numbers of distinct permutations 7r, which is 
given by a corresponding binomial coefficient. □ 



How can one use this result in the tomography entan- 
glement test of a PDC source? First each party measures 
along all three different polarization axes. Next one ei- 
ther actively post-processes the double click events or 
just passively computes the corresponding rates and/or 
probabilities of the full operators. Afterwards both par- 
ties can safely use the single photon assumption, or more 
precisely, the set of perfect single photon target operators 
{Ti^p} to compute the corresponding two-qubit state pab 
(single photon subspace on each side) via their preferred 
reconstruction technique. In case that this reconstructed 
state is entangled one can be assured that the observed 
data still verify entanglement if both parties believe in 
the more realistic measurement description {i 7 !^}. 

Next let us focus on the imperfections of the photo- 
detectors. Real photo-detectors register only some por- 
tion of the incoming photons, a significant part is not de- 
tected. If both detectors in the setup of Fig. [2] have the 
same inefficiency, this inefficiency can be modeled by an 
additional beam splitter in front of the perfect measure- 
ment device (3(|, hence if one combines the beam split- 
ter map (completely positive) with the already proven 
positive squashing map from the perfect case then one 
directly extends the validity of the positive squashing 
property to an inefficient measurement description. The 
same idea applies to dark counts, which can be modeled 
as a particular post-processing scheme on the classical 
outcomes (23|, and to misalignment errors, that arc de- 
scribed by a fixed depolarizing channel onto each photon 
separately. Even the extension to a multi-mode descrip- 
tion is possible if one employs the model from Ref. [37J ■ 

Concerning real experiments, one should note that 
double-clicks in a spatial mode can arise from different 
physical mechanisms. First, it can happen that due to 
the higher orders in the PDC process more than the de- 
sired number of photons arc generated and injected into 
the setup. Second, dark counts may lead to double click 
events. Finally, double-click events may occur in special 
setups for the generation of certain multipartite states, 
this case is, however, not important for our discussion 

Then, it is worth mentioning that the post process- 
ing used in the above scheme is usually not applied in 
real experiments: double click events are typically just 
thrown away. In practice, however, the amount of these 
undesired events is quite small: For instance, in the 4- 
photon experiment of Ref. [l9| the number the coinci- 
dences where a double click occurs in one mode while in 
the other three modes there is one click, is around 0.77 
% of the (desired) events, where in each mode there is 
exactly one click [59|]. It should be noted, however, that 
in experiments with more and more photons, these rates 
can be higher [38| , so that the penalty effect of the post- 
processing scheme becomes higher. 

Additionally we comment on two points: As one can 
prove, the corresponding squashing map is completely 
positive on the single and two photon subspace [l?] ]. 
Hence one only observes a violation of positivity of the 
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corresponding target density operator if the local multi- 
photon contributions are very large in comparison to the 
single and double photon part (and even then only for 
very particular entangled states); consequently, to ob- 
serve such a non-positive target operator in a real PDC 
experiment is very unlikely . 

As a last point we should make it clear that one can- 
not always apply Theorem IV. 11 Especially in multipar- 
tite experiments, it happens that one does not even want 
to obtain full tomography onto the multipartite target 
space but instead tries to measure an entanglement wit- 
ness with the least number of different global measure- 
ment settings. This may require more than three differ- 
ent settings on each photon. For instance, in the six- 
photon experiment of Ref. (39j an entanglement witness 
was measured which required seven measurements set- 
tings of the type Mi <g> Mi ® . . . <£> Mi , which is a significant 
advantage compared with the 3 6 = 729 settings required 
for state tomography. However, on each photon seven po- 
larization measurements have been made and the target 
observables are tomographically overcomplete. In such 
cases this theorem does not apply, because the linear de- 
pendencies imposed by the target operators are not satis- 
fied by the full observable set, cf. Eq. ([3]), hence the local 
squashing operation does not exist — in fact the map can- 
not even be linear. Here one might attempt to proceed 
with a global, separable squashing operation, cf. com- 
ment after Prop. IIII.ll 



VI. POSITIVE SQUASHING AND 
ENTANGLEMENT QUANTIFICATION 

In this section we argue that a local squashing opera- 
tion, even if it is not completely positive, can in principle 
not only provide qualitative indications about the pres- 
ence of entanglement, as was proved in Proposition lIII.ll 
but also quantitative ones. 

Let us start by recalling the notion of entanglement 
measure. An entanglement measure is a function from 
density operators to (positive) real numbers, that cap- 
tures quantitatively some property of entangled states. 
There are many entanglement measures in the litera- 
ture [40| ; some of them have an operational character, 
while some others focus on structural properties of en- 
tangled states, for example, the fact that, by definition, 
they do not admit a separable decomposition. Even if 
some entanglement measures do not have a direct opera- 
tional interpretation, they arc nonetheless useful because 
they may provide upper and lower bounds to operational 
measures or other interesting quantitative properties of 
entanglement. Furthermore, any entanglement measure 
can be considered as a benchmark for the quality of an 
experiment designed to create "highly entangled" states 
and to display a good global control on more than one 
subsystem at a time. This is due to two facts. The first 
is that, although different entanglement measures do not 
correspond to the same ordering of states from "unen- 



tangled" to "maximally entangled", there is typically a 
correlation: a state that is highly entangled with respect 
to one measure is, in most cases of interest, highly en- 
tangled with respect to another one. The second fact is 
that in an axiomatic approach to entanglement measures, 
it is typically asked that entanglement, as quantified by 
some entanglement measure, does not increase under the 
restricted class of Local Operations and Classical Oper- 
ations (LOCC). Indeed, entanglement cannot be created 
by LOCC operations alone, and it is natural to require 
that any entanglement measure does not increase under 
this set of operations. In this way, on one hand, entan- 
glement is elevated to a resource that by LOCC can only 
be manipulated and not augmented, and on the other 
hand, entanglement measures satisfying such an LOCC 
monotonicity are a fair benchmark for the ability of the 
experimenters to jointly manipulate many subsystems. 

Let us be more precise about the LOCC monotonicity 
of entanglement measures, focusing on the bipartite case. 
We say that E is an LOCC monotone if 

E{ P ab) > £(A LO cc[pab]), (15) 

where Alocc is an LOCC transformation. In particular, 
local completely positive trace-preserving (CPTP) maps 
belong to the class of LOCC operations, so that E is 
monotone with respect to CPTP local operations: 

E{ P ab) > E((A A ® A b )[pab])- (16) 

Thus, if the squashing is realized by local CPTP maps, 
the entanglement of the reconstructed squashed state 
(Aa <S> Ab)[pab] is a lower bound for the entanglement 
of the physical state actually prepared. The point here 
is that one can generalize Eq. (fT6)) to the case of posi- 
tive but not completely positive maps, at least for the 
entanglement measure called negativity [4l], , which is 
one of the few entanglement measures that can be easily 
computed. 

The negativity of a bipartite state pab is defined as 

N(p) = Mstzl, (17) 

where p AB = (J 1 ® id)[pAB] denotes the partial transpose 
of the original density operator. Here, T is the trans- 
position, which is a positive but not completely positive 
trace-preserving map, while "id" stands for the identity 
map, and ||A||i = tr(vAtA) is the trace norm on oper- 
ators. The value of the negativity is independent of the 
party we choose to apply transposition to, and quantifies 
the degree of violation of the partial transposition sepa- 
rability criterion [ill |44|, Indeed, it corresponds to the 
sum of the absolute values of the negative eigenvalues of 
Pab- 

In the Appendix we prove the following inequality 
N(pab) > -r 1 — T 77 [N((A A ® A b )[pab]) 

||Aa®a b ||F\ 

llA A ®A B ||f -n 
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with Aa = To Aa o T being (completely) posi- 
tive if and only if Aa is (completely) positive, and 
with a norm on linear maps defined as = 
ma^ ):W)=1 ||Q[|V>>(V>|]||i (cf- Ref. [H). We stress 

that || Aa ® Ab || is a measure of the joint violation of 
complete positivity of Aa and Ab- Indeed, if both maps 
Aa and Ab are trace-preserving and completely positive 
then one obtains ||Aa <S> AbU^ = 1 and one recovers the 
inequality given by Eq. (Tl6l) . 

Let us remark that N((Aa <£> Ab)[pab]) is the nega- 
tivity, as defined by Eq. (fT?)) , of the Hermitian operator 
(Aa® Ab) [pab] • If the local squashing operations Aa and 
Ab are not completely positive, then the latter need not 
be a physical state because of negative eigenvalues even 
before the partial transposition. The correcting terms in 
Eq. (p~8|) . with respect to Eq. (fl6|) . in particular the pres- 
ence of || Aa® Ab||^, take care of this possibility, making 
inequality (fT5|) hold. 

As shown in the Appendix, a different and possibly 
weaker lower bound on the negativity is given by 

N(pab) > .. .,\ jj (iV((A A ® A b )[pab]) 

||Aa||o||A b ||o V / iq x 

||A A ||o||A B ||o — 1 ^ 1 ' 

Here ||f2||<> = ||f2<8>id||i is the diamond norm 46], with 
the identity map that can be considered as acting on the 
same input space as fi, and ||f2||i = sup|| X || 1<1 ||f2[X]||i 
the trace norm for maps. 

We further remark that, in the case of a positive but 
not completely positive squashing operation, it might be 
possible to obtain lower bounds for the entanglement of 
the original state also for other entanglement measures. 
Although we are unable to provide further explicit exam- 
ples at this time, we observe that this might be true for 
the relative entropy of entanglement [13, |48| . The latter 
is defined for a state pab as 

Er{p ab) = min5(pAB|kAB), (20) 

where the minimum runs over all separable states and 
•SQcabIIo-ab) = tr[pAB(log 2 pas, - log 2 (7AB)] is the rela- 
tive entropy. Monotonicity of this measure under CPTP 
LOCC operations can be easily checked as follows: 

E r(p ab) = minS(pABljoAB) 

> min 1 S , (A L occ[PAB]||A L occ[o-AB]) 

(21) 

> min5 , (A L occ[/3AB]||TAB) 
tab 

= ^(AloccIpab])- 

In the first inequality one uses monotonicity of the rela- 
tive entropy under CPTP maps; for the second inequality 
one employs the fact that a CPTP LOCC map transforms 
a separable state into another separable state. Now, a lo- 
cal map Aa <8> Ab also transforms a separable state into a 



separable state as soon as the maps Aa and Aa are pos- 
itive and trace-preserving — this was the key fact used in 
Proposition If 11. 1 1 If monotonicity of the relative entropy 
holds under some local map Aa <8> Ab , even if Aa and Ab 
are just positive but not completely positive, then the 
inequality Er(pab) > Er((Aa <8> Ab)[pab]) still remains 
true. This possibility is left open for example by the fact 
that the requirement often used to prove monotonicity 
of the relative entropy is not c omp lete positivity, but the 
weaker request of 2-positivity [4!| (together with a trace 
preservation condition). A given map Q is 2-positive if 

where A, B, C and D are matrices themselves. Hence, if 
both maps Aa and Ab are positive and trace-preserving, 
and the combined local map Aa <£> Ab is 2-positive, then 
the inequality E r (pab) > E R ((A A ® A b )[pab]) still 
holds. 

In conclusion, a positive squashing operation does not 
only provide qualitative statements about entanglement, 
but potentially also quantitative ones. Open problems 
regard the application of the derived bounds on the neg- 
ativity to specific cases, and the analysis of other en- 
tanglement measures. Let us remark that in case of the 
negativity a detailed analysis of lower bounds on the en- 
tanglement essentially deals with the issue of separating 
the negativity due to the application of the local squash- 
ing maps from the negativity due to partial transposition. 
As the bounds are conservative, only states that are suffi- 
ciently entangled may result in a non-trivial lower bound. 
Indeed, it is clear that if AT((Aa<£> Ab)[/9ab]) = — this is 
the case for a separable pab and positive Aa and Ab — 
and ||A A ® A B ||f > 1 or || A A ||o|| A B || > 1, respectively, 
then the right-hand sides of Eq. p8[) and Eq. (|19[) are ac- 
tually negative. It is worth remarking that in the deriva- 
tion of the bounds for the negativity we have not made 
use of the positivity of the squashing operations. This 
indicates that if one considers local squashing operations 
with the aim of entanglement verification and quantifica- 
tion, then one may hope to be able to further relax the 
requirements on the corresponding maps, not only going 
beyond complete positivity, but also beyond positivity if 
adequate care is taken. 

VII. CONCLUSIONS 

Entanglement verification typically assumes that one 
knows the underlying measurement operators so that 
each classical outcome gets an accurate quantum me- 
chanical interpretation. We have addressed the question 
under what conditions an affirmative entanglement state- 
ment remains valid if a simplified description of the mea- 
surement apparatus is used. This situation can occur 
if the actual measurement observables are different from 
the ones used in the verification analysis, simply because 
of imperfections or wrong calibration. However it can 
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even happen on purpose: Indeed one can try. despite 
being aware of certain differences, to explain the data 
via an oversimplified model, e.g., a very low-dimensional 
description, that eases then the task of applying an en- 
tanglement criteria. Such a case occurs for example for 
an active polarization measurement with threshold de- 
tectors to analyze the entanglement from a PDC source. 
Here one may choose a single photon description only, 
although one knows that certain multi-photon states can 
also trigger events that are indistinguishable from a single 
photon case, because then one easily obtains "tomogra- 
phy" by using three different measurement settings and 
directly checks for entanglement on the reconstructed 
two-qubit state. 

Summarizing, a positive entanglement statement re- 
mains valid if the two operator sets can be related by 
a positive (not necessarily completely positive) map. In 
case that the reconstruction operators provide complete 
tomography such a positive maps exists if and only if 
all measurement results from the refined, actual mea- 
surement devices are compatible with the assumed mea- 
surement description of the device. We have shown that 
the aforementioned polarization measurement, measured 
along all three different polarization axes, constitutes a 
physical relevant example of such a connection that is 
positive but not completely positive. This result shows 
that most of the performed tomography entanglement 
tests for a PDC source are indeed error-free if one incor- 
porates a penalty for double clicks during the state recon- 
struction. This verifies entanglement for a more realis- 
tic model, with imperfections and multi-photon contribu- 
tions, of the measurement used. Finally, we argued that 
it might be possible to obtain not only a positive qualita- 
tive statement about the presence of entanglement, but 
also a quantitative one, even in cases where the squash- 
ing map is not completely positive and standard results 
about monotonicity of entanglement measures can not 
directly be used. 
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APPENDIX 

For a Hcrmitian operator pab normalized to satisfy 
tr(pAB) = 1, we define negativity as 



where p^ B — (T <E> id) [pab], and T is the transposition. 
Negativity corresponds to the sum of the negative eigen- 
values of Pab- 

Any Hcrmiticity preserving map A can be written as 
A[X] = ^CiKiXKj, c, e R. Then To A = A o T, 
with A : X h->- CiK*XKj, i.e., A = T o A o T. If A 
is (completely) positive, that is if Cj > for all i, then 
A is (completely) positive. It also holds that A is trace- 
preserving if and only if A is trace-preserving. 

For any map f2 we define the norm ||f2||^ = 

max |^):<V|V>=ill fi ^^^l 1 E^' Moreovcr ' we observe 
that the trace norm of a Hcrmitian operator X can be 
expressed as ||A||i = max_j <M< j tr(MA). For any 
-1 < M < a, 

|^|Ot[M]|^| = |tr(Ot[M]|^)^|)| 

= |tr(Mfi[h/>)(V|])| (24) 

<l|fi||f- 

Therefore, if -1 < M < 1, then -1 < < 1. 

Thus, assuming that Aa and Ab are trace- preserving — 
so that tr((A A ® A b )[pab]) = 1: 

A((A A <8> A b )[pab]) 

= ||((roA A )® A b )[pab]||i - 1 
2 

_ ||(A A <8> Ab)[p ab ]IIi -1 



-{ max tr(Af(A A ® A B )[p AB ]) -1 

2 L -1<M<1 

-\ max tr((A A ® A B ) t [M]p AB ) - 1 

2 I -1<A/<1 



-{||Aa® A B ||f max tr 

2 t -1<M<1 



(Aa® A B )t[M] r 
||A A ®A B ||f PAB 



< 



i{||A A <g>A B ||f max tx(M'p r AB )-l\ 

{||AA®AB||f||p AB ||-l} 

\A A ® A B ||f A(p AB ) + - 



A a ® A 



Bill 



(25) 



Solving for A (pab) we finally find 
1 



N(pab) > 



|A A ® A B ||f 



A((A A ® A b )[pab]) 



\A A ® A B ||f - 1 



(26) 



N(pab) 



IpLIIi-i 



(23) 



For a Hcrmiticity preserving map one easily checks 
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that o 9||f < ||ft||f ||6||f . In our case 



||A A ® A B ||f < ||AA®id Bout ||f ||id A „ 



A 



B||f 



< 



\Aa ® ide 

|a a ||o||Ab 



Ifl 



idA h 



A B ||f 



(27) 



By idB out and idA in we have denoted the identity map 
on the output space of Ab and on the input space of Aa, 
respectively. The equality in the second line is due to the 
fact that ||ToftoT||f=||ft||f. The diamond norm [H 
is defined as ||ft|| = ||ft<£> id||i, with the identity map 
that can be taken as acting on the same input space as 



ft, and with ||ftj|i = sup|| X jj i<:L |Jft[A]||i the trace norm 
for maps. The last inequality is due to the fact that 
||ft®id||f < ||ft <g> idjji < ||ft||o, for id acting on an 
arbitrary dimension [46J|. Thus, we finally obtain the 
lower bound 



^(pab) > 



A 



A\\o\ As \\o 



iV((A A g> A b )[/9ab]) 



I A 



A||o 



IA, 



- 1 



(28) 
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